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ABSTRACT 

Flow  visualization  of  low  speed  flow  has  been  obtained 
in  this  thesis  by  the  injection  of  a  thin  stream  of  gas. 
This  injected  gas  has  a  refractive  index  which  differs  from 
that  of  the  uniform  flow  in  the  wind  tunnel.   A  frozen  fringe 
hologram  of  the  introduction  of  this  gas  indicates  regions 
of  fringe  shift.   A  pathline  is  produced  by  connecting  points 
of  maximum  fringe  shift.   Pathlines  generated  in  such  a  man- 
ner around  a  cylinder  agree  with  the  stream  functions  of 
potential  flow  theory. 

A  rapid  means  has  been  developed  of  data  reduction  of 
fringe  information  from  a  photograph  of  a  reconstructed  holo- 
graphic interf erogram.   Data  points  are  gathered  using  a 
digitizer  and  a  mini-computer.   The  array  of  data  points  so 
obtained  can  be  examined  and  corrected  if  necessary  and  fin- 
ally expanded  by  means  of  Lagrangian  interpolation.   The 
expanded  array  is  transferred  to  a  fortran  program  which  ac- 
complishes the  integral  inversion  and  subsequently  prints 
out  the  density  field.   By  using  this  method  of  data  reduc- 
tion and  inversion,  a  large  time  savings  is  realized. 
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I.   INTRODUCTION 

This  thesis  consists  of  two  parts.   The  first  part  is  an 
investigation  of  low  speed  flow  visualization  by  means  of 
holographic  interf erometry .   The  lack  of  significant  density 
changes  over  a  profile  at  low  speeds  has  made  interf erometric 
investigations  of  low  speed  flow  difficult.   In  this  study, 
flow  visualization  at  low  speeds  has  been  obtained  by  inject- 
ing a  thin  stream  of  gas  in  the  wind  tunnel.   This  gas  had  a 
refractive  index  which  varied  from  that  of  the  uniform  flow 
in  the  wind  tunnel.   A  frozen  fringe  hologram  of  the  flow 
field  with  the  introduced  gas  indicated  regions  of  fringe 
shift.   The  pathline  can  be  found  by  connecting  points  of 
maximum  fringe  shift.   Pathlines  generated  in  such  a  manner 
were  compared  with  potential  flow  theory  by  superimposing  a 
family  of  stream  functions  over  the  experimental  pathline. 
This  comparison  was  done  with  a  cylinder  profile  in  a  uniform 
steady  flow. 

The  second  part  of  this  thesis  concerns  itself  with  a 
rapid  means  of  data  reduction  of  the  large  amount  of  fringe 
data  that  is  obtained  in  holographic  interf erometry .   A  dig- 
itizer and  a  mini-computer  were  used  to  effect  this  rapid 
means  of  data  reduction.   The  number  of  data  points  is  ex- 
panded in  the  mini-computer  by  using  Lagrangian  interpolation 
The  whole  array  of  points  are  also  automatically  transferred 
to  a  fortran  data  reduction  program.   This  existing  fortran 
program  then  inverts  the  Integral  equations  to  determine  the 


density  field.   In  the  past,  fringe  data  had  been  collected 
from  a  hologram  by  the  use  of  an  optical  comparator  to  deter- 
mine the  amount  of  fringe  shift  along  each  plane.   The  data 
was  typed  onto  IBM  cards  and  fed  into  a  fortran  program  for 
reduction.   The  present  development  of  the  use  of  a  digitizer 
and  the  ability  to  automatically  feed  data  into  a  fortran 
program  greatly  reduces  the  time  to  obtain  a  density  field 
from  a  hologram. 
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II.   FLOW  VISUALIZATION 

A.   HOLOGRAPHY  AND  FLOW  VISUALIZATION 

Dennis  Gabor  in  1948  (1)  was  the  first  to  invent  a  meth- 
od of  wavefront  reconstruction  called  holography.   With  the 
development  in  the  early  1960's  of  a  coherent  source  of  ra- 
diation, the  laser,  the  field  of  holography  rapidly  developed. 

Holography  is  the  process  of  storing  wavefronts  in  a 
photosensitive  material  called  a  hologram.   At  a  later  time 
these  wavefronts  may  be  reconstructed  from  the  hologram  to 
form  an  image  which  has  the  same  depth  and  parallax  charac- 
teristics as  the  real  object.   This  is  accomplished  by  re- 
cording the  interference  pattern  between  the  light  diffracted 
by  the  object  and  a  reference  beam  of  light.   An  optical 
arrangement  is  set  up  to  split  the  coherent  light  source  into 
a  reference  beam  and  a  scene  or  object  beam.   Within  certain 
restrictions  the  paths  of  these  two  beams  are  arbitrary. 
The  photosensitive  plate  is  placed  in  such  a  position  that 
the  reference  beam  and  the  diffracted  light  from  the  object 
fall  upon  it.   As  stated  in  Ref.  1,  each  point  on  the  holo- 
gram contains  the  full  information  of  the  object  as  seen 
from  that  point.   Full  information  constitutes  intensity  and 
phase  information.   The  phase  information  is  derived  from 
the  interference  pattern.   When  the  object  is  reconstructed, 
each  of  these  points  superimposes  its  own  complete  picture, 
which  has  depth  and  parallax.   The  recording  of  phase  infor- 
mation (with  the  associated  dc   h  and  parallax)  i   'hat 
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distinguishes  holography  from  photography  which  only  records 
intensities . 

A  frozen  fringe  hologram  is  one  which  has  been  doubly 
exposed.   One  exposure  records  the  reference  situation.   The 
second  exposure  records  whatever  perturbation  has  been  intro- 
duced.  In  the  present  low  speed  flow  visualization  experi- 
ment the  reference  condition  of  exposure  is  the  steady  flow 
of  uniform  gas  over  the  two-dimensional  profile.   After  the 
reference  exposure  an  optical  component  is  moved  slightly  to 
provide  fringe  lines.   The  flow  field  is  perturbed  by  the 
introduction  of  a  thin  stream  of  gas  into  the  uniform  steady 
flow  ahead  of  the  profile.   A  second  exposure  is  taken  on 
the  same  hologram  plate  of  the  perturbed  flow.   In  the  re- 
construction of  the  resulting  hologram,  the  two  recorded 
flow  fields  interfere  with  each  other  and  show  a  deviation 
in  the  superimposed  finite  fringes.   Since  both  exposures 
utilize  the  same  optical  components,  the  fringe  lines  are 
lines  of  constant  differences  in  optical  path  lengths.   Any 
imperfection  in  the  optical  quality  of  the  apparatus  is  thus 
eliminated. 

When  a  gas  of  a  different  density  (optical  path  length) 
is  carefully  introduced  into  the  flow  field,  it  is  swept 
through  the  wind  tunnel  and  follows  a  pathline.  Holography 
is  extremely  sensitive  to  variations  in  optical  path  length 
and  by  using  holographic  techniques  a  line  of  fringe  shifts 
can  be  seen  in  the  region  of  the  injected  gas.  By  carefully 
adjusting  the  amount  of  gas  introduced  in  relation  to  the 
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wind  tunnel  velocity,  a  small  well  defined  line  of  fringe 
shifts  can  be  produced.   A  pathline  can  be  created  by  con- 
necting points  of  maximum  fringe  shift. 

B.   EXPERIMENTAL  APPARATUS 

The  arrangement  of  the  experiment  is  shown  in  figure  1(a). 
This  type  of  set-up  is  known  as  a  modified  Z-Type  Schlieren. 
A  Korad  Q~switched  laser  (ruby  wave  length  6943  A)  is  mounted 
on  an  optical  rail.   A  small  five-milliwatt  helium-neon  con- 
tinuous gas  laser  is  placed  adjacent  to  the  rail  to  provide 
alignment  of  the  optics.   The  beam  is  expanded  and  projected 
onto  a  parabolic  mirror.   This  mirror  collimates  the  beam  to 
facilitate  easier  reconstruction.   The  reference  beam  travels 
through  the  beam  splitter  and  is  reflected  off  two  mirrors 
before  illuminating  the  holographic  plate.   The  object  beam 
is  reflected  off  the  beam  splitter  and  then  projected  through 
the  wind  tunnel.   Just  before  the  tunnel  the  object  beam 
passes  through  a  diffusing  plate  and  a  grid.   The  actual  ex- 
posure time  with  the  Q-switched  laser  was  on  the  order  of 
twenty  nano-seconds.   This  stopped  all  movement  during  the 
exposure,  but  did  not  prevent  movement  between  exposures. 
The  tunnel  was  mounted  separately  from  all  other  components 
to  minimize  unwanted  movement  caused  by  vibration.   Alignment 
of  the  optics  was  done  by  projecting  the  small  laser  beam 
through  the  ruby  laser  head  and  out  into  the  optical  system. 
The  initial  alignment  of  the  small  laser  with  the  Q-switched 
laser  was  accomplished  with  burn  holes  in  photographic  paper. 
This  technique  is  explained  in  Appendix  A  of  Ref.  2.   Fringe 
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lines  were  formed  by  moving  the  beam  splitter  .003  inches  be- 
tween exposures.   The  beam  splitter  was  mounted  on  a  preci- 
sion X-Y  translating  stage. 

The  wind  tunnel  is  3.5  inches  by  3.5  inches  and  construct- 
ed out  of  plexiglass  (not  of  optical  quality).   The  speed  of 
the  tunnel  can  be  varied  from  zero  to  twenty-two  feet  per 
second.   A  section  of  honeycomb  was  placed  at  the  entrance  to 
serve  as  a  flow  straightener  and  as  a  support  for  the  tube 
which  injects  the  gas. 

The  tube  used  for  injecting  the  gas  had  an  outside  dia- 
meter of  3/32  inches.   It  was  connected  by  tubing  to  a  flow 
meter.   This  meter  allowed  visual  adjustment  of  the  injection 
gas  flow  rate. 

Two  profiles  were  used,  both  constructed  from  plexiglass. 
Pictures  of  the  profiles  and  their  supports  are  shown  in  Fig. 
1(e).   A  cylinder  was  used  having  a  length  of  three  inches 
and  a  diameter  of  one  inch.   The  second  profile  was  a  wing 
having  a  length  of  three  inches,  a  chord  of  1.5  inches  and  a 
thickness  of  3/16  inches.   Both  profiles  were  supported  by 
two  upright  pieces  of  plexiglass.   A  grid  was  lightly  inked 
onto  the  supports  to  provide  a  reference  during  reconstruc- 
tion.  The  ends  of  the  profiles  were  painted  black  for  better 
definition  of  thoir  outline.   There  was  a  snug  fit  between 
the  pins,  which  extended  from  each  end  of  the  profiles,  into 
the  two  supports.   This  fit  allowed  the  wing  to  be  set  at  a 
desired  angle  of  attack  and  remain  there.   When  the  profile 
was  adjusted  between  the  supports  then  the  assembly  was  slid 
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into  the  tunnel.   The  honeycomb  was  placed  at  the  entrance 
and  the  injection  tube  inserted  through  the  honeycomb. 

Reconstruction  was  done  with  a  fifteen-milliwatt  helium- 
neon  laser.   The  beam  was  expanded  from  .04  inches  to  .25 
inches  by  the  use  of  two  mirrors.   A  sketch  of  the  recon- 
struction is  shown  in  Fig.  1(d).   A  shutter  was  placed  be- 
tween the  second  mirror  and  the  hologram  plate  holder.   In 
order  to  photograph  the  real  image  the  camera  lens  was  re- 
moved.  A  dark  cloth  was  used  to  cover  the  laser  and  prevent 
any  extraneous  light  from  overexposing  the  film. 

Other  equipment  used  in  this  experiment : 

1.  power  supply  and  circulator  for  Q-switched  laser 

2.  eye  protectors 

3.  Agfa  10E75  holographic  emulsion  plate  and  Polar- 
oid Type  55  film 

4.  auto  collimator  for  internal  Q-switched  laser 
alignment . 

C.   EXECUTION  AND  RECONSTRUCTION 

Before  taking  a  frozen  fringe  hologram,  the  alignment  of 
the  optical  system  was  verified  using  the  small  5-milliwatt 
laser.   While  this  was  being  done  the  power  supply  for  the 
Q-switched  laser  was  turned  on.   The  recommended  warm-up 
time  for  the  power  supply  is  five  minutes.   The  circulator 
was  also  turned  on  to  stabilize  the  laser  head  and  the  etalon 
at  a  low  temperature.   A  small  amount  of  air  was  blown  onto 
the  optical  surfaces  to  remove  dust.   The  mirror  directing 
the  alignment  beam  through  the  ruby  was  removed  and  a  2.5mm 
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transverse  mode  selector  pin  hole  was  placed  in  its  holder. 
The  Agfa  10E75  emulsion  plate  was  placed  in  the  plate  holder 
which  was  under  a  dark  cloth.   After  appropriate  warnings  to 
other  personnel  in  the  area  and  donning  safety  glasses,  the 
laser  was  charged.   The  wind  tunnel  was  set  at  the  desired 
speed  and  allowed  to  stabilize.   The  black  cloth  covering 
the  holder  was  removed,  the  laser  fired,  and  the  black  cloth 
replaced.   Then  the  beam  splitter  was  moved  .003  inches. 
The  injection  gas  was  turned  on  and  the  flow  regulated.   The 
capacitor  bank  for  the  laser  was  again  charged  and,  after 
removing  the  cloth,  fired.   The  laser,  wind  tunnel  and  the 
injection  gas  were  turned  off  while  developing  the  emulsion 
plate. 

Emulsion  plate  development  was  standard: 

1.  Kodak  D-19  developer  for  five  minutes 

2.  acetic  acid  stop  bath  for  thirty  seconds 

3.  fixer  for  five  minutes. 

The  plate  was  washed  in  water  followed  by  a  rinse  in  Kodak 
Photoflo  wetting  agent.   A  hot  air  dryer  was  used  to  dry  the 
plate.   Occasionally  the  plate  was  overexposed  and  a  bleach 
was  used.   Whenever  this  was  done  an  index  matching  solution 
of  xylene  was  also  used  to  reduce  the  noise  level  of  the 
hologram.   Further  information  on  bleaching  and  index  match- 
ing is  available  in  section  10.8.5  of  Ref.  1. 

Reconstruction  was  done  with  a  fifteen-milliwatt  Stabi- 
lite  laser  using  the  setup  explained  in  section  B   and 
illustrated  in  Fig.  1(d).   The  virtual  Image  was  quickly 
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reconstructed  to  assess  the  quality  of  the  hologram.   This 
was  easily  done  by  using  the  collimating  lens  of  the  Stabi- 
lite  laser  and  looking  through  the  hologram.   If  the  quality 
of  the  plate  warranted  a  photographic  record  then  the  plate 
was  reversed  in  its  holder.   The  collimating  lens  was  re- 
moved from  the  laser  and  the  beam  was  slightly  expanded  to 
reduce  speckle.   The  virtual  image  was  projected  onto  the 
frosted  glass  at  the  back  of  a  lensless  camera.   For  a  dis- 
cussion of  real  and  virtual  images  the  reader  is  directed  to 
Ref.  1.   An  f-16  aperature  was  used  on  the  shuttering  device. 
Black  cloth  covered  the  camera,  holder,  shutter  and  the  la- 
ser.  The  plate  holder  was  adjusted  to  present  the  correct 
picture  on  the  ground  glass  of  the  camera.   A  photographic 
record  was  made  using  Polaroid  Type  55  positive/negative 
film. 

D.   COMPARISON  OF  RESULTS  WITH  POTENTIAL  FLOW  THEORY 

Two-dimensional  potential  flow  theory  was  used  to  give  a 
comparison  of  experimental  pathlines  obtained  by  connecting 
points  of  maximum  fringe  shift.   A  discussion  of  this  theory 
is  presented  in  Ref.  4. 

The  program  listed  in  Appendix  A  calculates  stream  func- 
tions of  uniform  flow  past  a  circular  cylinder  with  circula- 
tion.  The  stream  function  i[»  is  given  as: 

,{,  =  Uy  -  lA4  +  ~  ln(-) 

r  J  r     2ir    va 


where    U  =  uniform  velocity  ;     r  =  vx2  +  y?; 

a  =  cylinder  radius;       t  =  circulation. 
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An  example  of  the  graphical  output  of  a  cylinder  with  circu- 
lation is  shown  in  Fig.  2(c).   In  the  experimentation  for 
this  thesis  there  was  no  circulation  around  the  cylinder. 

When  the  photograph  of  the  cylinder  in  Fig.  2(a)  is  super- 
imposed over  the  stream  functions  calculated  from  the  program 
in  Appendix  A  then  Fig.  2(d)  is  the  result.   Notice  the  close 
agreement  between  the  stream  function  labeled  0  and  the 
points  of  maximum  fringe  shift  over  the  profile.   Downstream 
of  the  cylinder  the  shift  in  the  stream  function  0  is  be- 
lieved due  to  the  heavier  carbon  dioxide  used  as  the  injected 
gas.   As  the  wind  tunnel  speed  is  increased  the  point  of  di- 
vergence of  the  injected  gas  with  the  stream  function  is 
shifted  further  downstream.   It  was  found  that  for  very  low 
tunnel  speeds  the  fall  of  the  carbon  dioxide  occurred  at 
earlier  points.   At  wind  tunnel  speeds  of  22  feet  per  second 
the  line  of  fringe  shifts  caused  by  injecting  carbon  dioxide 
agreed  favorably  with  a  stream  function  over  a  cylinder. 

E.   DISCUSSION  OF  FLOW  VISUALIZATION 

Three  types  of  injection  gases  were  investigated.   Helium, 
argon,  and  carbon  dioxide  were  introduced  into  the  flow  field. 
Helium  and  carbon  dioxide  provided  an  acceptable  fringe  shift. 
The  most  advantageous  gas  was  carbon  dioxide  since  its  high 
molecular  weight  prevented  the  injected  stream  from  rapidly 
diffusing  in  the  tunnel.   Fig.  2(a)  shows  that  this  fringe 
shift  is  small  and  well  defined  over  the  profile. 
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To  accurately  define  the  pathline  the  injected  gas  must 
be  swept  through  the  tunnel.   If  the  flow  rate  out  of  the 
injection  tube  is  too  high,  then  the  subsequent  path  of  the 
gas  does  not  follow  a  pathline.   To  ensure  pathline  accuracy 
the  velocity  of  injection  gas  is  matched  as  closely  as  pos- 
sible to  the  velocity  of  the  tunnel.   Both  the  tunnel  veloc- 
ity and  the  injected  gas  velocity  were  22  feet  per  second. 
Velocity  matching  required  calibration  of  both  the  wind  tun- 
nel and  the  injection  flow  metering  device.   Once  the  veloc- 
ities are  matched  then  the  amount  of  gas  injected  is  dependent 
on  the  injection  tube  diameter.   At  low  speeds  a  tube  which 
had  a  diameter  5%  of  the  height  of  the  tunnel  was  able  to 
define  a  fringe  shift  throughout  the  hologram.   At  higher 
tunnel  speeds  a  smaller  tube  can  be  used  since  the  gas  flow 
rate  is  increased  and  more  gas  is  introduced  to  cause  a 
shift.   Fig.  2(b)  is  an  example  where  the  flow  velocities 
were  not  matched.   Notice  in  the  figure  how  the  injected  gas 
is  blown  onto  the  lower  surface  of  the  wing.   The  Kutta  Con- 
dition is  easily  recognized  in  the  circled  area  at  the 
trailing  edge. 

F.   CONCLUSION  DERIVED  FROM  FLOW  VISUALIZATION 

The  sensitivity  of  holography  to  slight  variations  in 
density  makes  this  flow  visualization  technique  ideal  for 
low  speed  flow  investigation.   Depending  on  the  nature  in 
which  you  inject  the  gas  one  can  emphasize  different,  fea  lures 
of  the  flow.   Although  it  was  not  done  for  this  thesis, 
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perhaps  vortices  can  be  visualized  by  injecting  more  gas. 
As  shown  by  Fig.  2(a),  when  the  velocity  of  the  injected  gas 
is  matched  to  that  of  the  wind  tunnel  the  resultant  line  of 
fringe  shift  makes  a  pathline  visible  around  the  profile. 
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III.   DATA  REDUCTION 

A.   RAPID  DATA  REDUCTION  TECHNIQUE 

Due  to  the  large  amount  of  data  generated  by  holographic 
interf erometry ,  a  rapid  means  of  data  reduction  is  desired 
in  the  determination  of  a  density  field  from  holography.   To 
date  the  emphasis  has  been  on  developing  a  better  fortran 
program  which  inverts  the  integral  equations  in  a  calculation 
of  a  density  field.   A  comparison  of  the  different  methods 
used  to  invert  the  integral  has  been  discussed  by  Sweeney 
[Ref.  6].   This  paper  addresses  the  problem  of  data  acquisi- 
tion and  transfer  of  these  data  to  a  fortran  inversion  pro- 
gram.  A  digitizer  is  used  in  determining  the  fringe  shift 
in  an  1 nterf erogram  obtained  by  reconstructing  the  frozen 
fringe  hologram.   A  mini-computer  processes  the  data  using 
Lagrangian  interpolation,  and  then  stores  it.   At  an  appro- 
priate time  the  data  are  recalled  from  the  storage  tape  cas- 
sette in  the  mini-computer  and  automatically  transferred  to 
a  large  computer  for  density  calculation.   In  the  previous 
method  of  data  reduction  the  fringe  data  had  been  collected 
by  the  use  of  an  optical  comparator  which  determined  the 
amount  of  fringe  shift  along  each  plane  of  a  hologram.   The 
fringe  data  were  then  typed  onto  IBM  cards  and  submitted  to 
a  fortran  program  in  a  large  computer.   The  present  develop- 
ment of  the  use  of  a  digitizer  and  the  ability  to  store  data 
and  then  automatically  transfer  it  to  a  fortran  program 
greatly  reduces  the  time  needed  to  obtain  a  density  field. 

1!) 


To  obtain  a  photograph  of  fringe  data,  a  doubly  exposed 
hologram  is  taken  of  the  aerodynamic  phenomenon.   The  first 
exposure  is  of  the  reference  condition  which  usually  is  a 
no-flow  situation.   Flow  is  then  introduced  into  the  area  in 
question  and  a  second  holographic  exposure  is  taken.   Finite 
fringe  lines  are  established  on  the  hologram  by  a  slight 
movement  between  exposures  of  an  optical  component.   The 
holographic  plate  is  developed  and  then  the  scene  is  recon- 
structed using  a  laser  source.   A  photographic  record  can  be 
taken  of  the  real  or  the  virtual  image.   For  a  further  ex- 
planation of  holographic  interf erometry  the  reader  is  direct- 
ed to  Ref.  5. 

The  actual  integral  inversion  in  the  fortran  program  can 
be  done  in  many  ways.   A  comparison  of  different  methods  has 
been  discussed  by  Sweeney  [Ref.  6].   Some  of  these  methods 
are: 

1.  Representation  of  the  integral  by  orthogonal  poly- 
nomials 

2.  Fourier  Expansions 

3.  Fourier  Transforms 

4.  Fast  Fourier  Transforms. 

The  overwhelming  amount  of  data  generated  by  holographic 
interf erometry  requires  a  rapid  method  of  data  acquisition, 
evaluation,  storage,  and  transfer. 
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B.   EQUIPMENT  DESCRIPTION  AND  OPERATION 

This  data  reduction  technique  uses  a  Hewlett  Packard  dig- 
itizer (HP9864A),  calculator  (HP9830A)  and  printer  (HP9866A). 
A  photograph  of  the  apparatus  is  shown  in  Fig.  3.   An  acous- 
tic coupler  was  also  used  to  allow  the  calculator  to  act  as 
an  on-line  terminal  with  the  Naval  Postgraduate  School  IBM 
360  computer.   The  power  supply  for  the  digitizer  was  mounted 
under  the  table  and  the  printer,  as  can  be  seen  in  Fig.  3, 
is  resting  on  top  of  the  calculator. 

Data  were  introduced  into  this  equipment  with  the  use  of 
the  cursor  attached  to  the  digitizer.   A  button  is  on  the 
cursor  which  defines  the  origin  on  the  magnetic  table.   The 
coordinates  of  single  points  are  entered  by  depressing  but- 
ton "S" .   Continuous  points  can  be  entered  by  depressing  "C" . 
This  last  feature  was  not  used  in  the  program. 

The  programmable  calculator  uses  Basic  language  and  has 
an  initial  memory  capability  of  1760  words.   Additional  mem- 
ory can  be  added  externally  or  internally  to  the  calculator. 
A  tape  cassette  with  32,000  word  storage  capacity  can  be 
used  with  the  calculator.   Program  or  data  can  be  stored  on 
or  retrieved  from  this  cassette.   The  entire  program  is  on 
the  first  portion  of  this  cassette.   The  remainder  of  the 
magnetic  tape  is  used  to  store  information  for  up  to  nine 
pictures,  each  picture  containing  twenty  fringe  lines  and 
twenty  planes.   By  employing  a  Terminal  Read-Only-Memory 
(Term  ROM)  the  calculator  becomes  an  on-line  terminal  for  an 
IBM  360  computer.   While  in  the  Term  ROM  mode  the  calculator 
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retains  its  basic  operating  ability.   It  is  possible,  al- 
though not  advisable,  for  an  operator  to  execute  a  terminal 
command  and  while  the  IBM  operation  is  being  performed,  use 
the  calculator  locally.   Operation  of  the  calculator  as  a 
terminal  was  in  accordance  with  the  normal  terminal  procedure. 

C .   PROCEDURE 

The  procedure  of  using  this  method  of  data  reduction  is 
extremely  straightforward.   A  listing  of  the  programs  to 
transfer  data  is  given  in  Appendix  B.   Each  of  the  six  pro- 
grams is  stored  in  a  file  on  a  tape  cassette.   Each  program 
is  referred  to  by  the  cassette  file  in  which  the  program  is 
stored.   The  arrangement  of  the  programs  were  designed  to 
allow  an  operator  to  collect  and  expand  fringe  data  for  up 
to  nine  pictures.   At  a  later  time,  when  the  terminal  mode 
for  the  IBM  360  is  available,  then  the  stored  data  are  sent 
to  the  large  computer  for  reduction.   A  flow  diagram  of 
these  operations  is  shown  in  Fig.  4.   A  description  of  the 
procedure  used  in  the  data  reduction  follows. 

A  photograph  of  a  frozen  fringe  hologram  was  taped  to  the 
table  of  the  digitizer.   If  the  hologram  did  not  have  a  well 
defined  grid  than  a  thin  transparent  grid  was  taped  over  the 
photograph  to  provide  a  reference.   The  calculator,  printer 
and  digitizer  were  turned  on.   The  data  reduction  tape  cas- 
sette was  placed  in  the  calculator  and  the  command  "LOAD" 
was  executed.   The  program  was  started  at  point   A   on  Fig. 
4  by  executing  the  command  "RUN".   Statement  20  requested 
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the  operator  to  enter  the  time,  day,  month  and  year.   These 
numbers  were  introduced  into  the  calculator  by  typing  them 
onto  the  display.   All  numbers  were  separated  with  commas. 
The  display  was  executed  and  the  program  (HP  file  0)  contin- 
ued to  run.   HP  file  0  printed  the  time,  day,  month,  year, 
and  the  program  organization.   At  statement  190,  HP  file  1 
was  automatically  retrieved  from  the  cassette  and  loaded  in- 
to the  calculator.   HP  file  1  asked  for  the  orientation  of 
the  grid.   The  origin  and  an  arbitrary  second  point  on  the 
horizontal  axis  were  entered  by  depressing  the  cursor  "0" 
(for  the  origin)  and  then  the  cursor  "S".   By  lifting  the 
cursor  J  inch  off  of  the  table  the  established  origin  was 
erased  and  the  orientation  procedure  had  to  be  repeated. 
Whenever  the  calculator  is  ready  to  accept  a  data  point  from 
the  digitizer  a  soft  beep  is  heard.   The  calculator  asks  for 
the  number  of  fringe  lines  to  be  used  from  the  picture.   Next 
the  number  of  planes  in  the  picture  is  entered.   The  program 
is  able  to  accept  up  to  20  fringe  lines  and  20  planes.   The 
equation  of  each  fringe  line  is  entered  by  entering  two  arbi- 
trary data  points  on  the  line  in  the  linear  region.   After 
all  lines  were  entered  the  program  prints  the  mean  and  stand- 
ard deviation  of  the  slopes  of  all  lines.   Also  the  differ- 
ences from  the  mean  of  each  slope  are  printed.   The  operator 
makes  an  evaluation  as  to  the  validity  of  his  fringe  lines 
by  comparing  the  slopes.   If  any  slope  varied  greatly  from 
the  mean  then  all  fringe  lines  can  be  reentered.   An  example 
of  the  printout  is  shown  in  Appendix  C.   If  the  fringe  lines 
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are  satisfactory  then  "CONT"  was  executed  and  HP  file  2  is 
loaded  into  the  calculator.   The  picture  number  and  the  num- 
ber of  data  points  in  each  plane  are  entered.   As  each  data 
point  is  entered  the  elevation  and  fringe  shift  are  printed. 
At  conclusion  of  each  plane  the  array  of  elevations  and 
shifts  are  stored  on  a  cassette  tape  file.   The  first  number 
in  the  array  is  the  number  of  data  points  in  that  plane. 
The  second  number  in  the  array  is  the  elevation  of  the  first 
data  point.   The  third  number  is  the  fringe  shift  for  the 
first  data  point.   Elevations  and  shift  alternate  through 
the  remainder  of  the  array. 

The  reference  file  is  filled  after  acquiring  all  the  data 
points  in  the  picture.   The  reference  file  must  be  compatible 
with  the  fortran  data  reduction  program.   A  program  based  on 
Fourier  Transform  methods,  [Ref.  3],  was  modified  to  accept 
data  input  parameters  in  a  specific  order.   After  the  refer- 
ence file  is  filled,  it  is  stored  in  the  cassette  under  the 
appropriate  picture  number.   The  parameters  in  the  reference 
file  are  shown  in  Appendix  B,  file  2  starting  with  statement 
430.   The  operator  was  at  point   B   in  Fig.  4.   Either  he 
can  gather  more  data  from  different  pictures  or  expand  the 
data  for  fortran  reduction. 

HP  file  4  is  loaded  for  data  expansion.   Expansion  is 
done  using  Lagrangian  interpolation.   All  of  the  data  points 
from  the  plane  to  be  reduced  are  printed  out.   The  program 
asks  for  the  data  point  numbers  which  will  be  the  basis  of 
the  array  of  expanded  data.   The  number  of  points  to  use  for 
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curve  fitting  should  be  less  than  seven.   If  a  greater  number 
of  points  are  used  then  the  computational  time  greatly  in- 
creases but  more  importantly  the  curve  fit  to  the  data  points 
is  erratic.   If  more  than  one  picture  needs  to  be  reduced 
then  the  same  point  numbers  are  automatically  used  for  data 
expansion  of  each  picture.   The  first  number  in  the  expander 
array  is  the  number  of  data  points  to  follow.   The  second 
number  is  the  number  of  the  plane  in  the  picture.   Before  the 
expansion  of  the  array  all  positions  were  initialized  to 
9.99999.   This  is  shown  in  File  FT03F001  of  Appendix  E  where 
the  actual  data  that  has  been  fed  in  are  indicated  and  those 
positions  in  the  array  that  are  not  expanded  data  are  9.99999. 
The  reason  for  initializing  the  array  is  for  convenience  of 
data  transmission.   The  absence  of  the  nines  following  the 
array  should  be  investigated  by  reading  out  the  stored  data. 
After  the  expanded  data  points  are  stored  in  the  cassette 
the  operator  loads  HP  file  3. 

File  3  is  the  program  which  automatically  feeds  the  data 
into  the  fortran  system.   Directions  for  its  use  are  printed 
by  executing  "LIST  10,110".   To  reduce  the  transmission  time 
the  expanded  data  is  sent  to  the  IBM  360  in  arrays  of  six 
numbers.   The  arrays  of  six  numbers  are  placed  on  a  storage 
disk  in  the  IBM  360  system  under  a  file  labeled  "FILE  FT03F001" 
The  fortran  reduction  program  calls  FILE  FT03F001  for  its 
data  input. 
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D.   DISCUSSION 

The  previously  described  technique  permits  extremely 
rapid  data  acquisition  and  transfer.   The  ability  to  rapidly 
acquire,  transfer,  and  store  fringe  data  in  large  quantities 
allows  the  operator  to  process  many  holograms.   To  determine 
a  density  field  around  an  asymmetric  object  in  a  wind  tunnel, 
as  many  holographic  views  as  possible  are  used.   Each  of 
these  holograms  is  taken  at  a  different  angle  and  provides 
fringe  data  as  seen  from  that  angle.   This  program  has  the 
ability  to  process  nine  holograms  of  400  data  points  each. 

An  improvement  of  this  data  reduction  method  can  be  made 
by  incorporating  a  plotter  peripheral  and  by  developing  a 
better  interpolating  routine.   At  present  the  expanded  data 
is  based  on  six  fringe  shift  data  points.   When  more  than 
six  points  are  used  the  expanded  array  becomes  erratic.   One 
would  like  to  utilize  as  many  of  the  available  20  data  points 
as  possible  on  each  plane.   A  new  interpolating  routine  could 
be  developed  to  accept  more  than  six  points  and  still  yield 
accurate  results.   The  expanded  array  and  the  original  data 
points  could  be  compared  on  the  plotter.   From  the  graphical 
output  of  the  plotter  an  improved  subset  of  data,  points  could 
be  chosen  to  expand  the  array.   By  incorporating  this  scheme 
not  only  the  best  order  polynomial  is  used  to  expand  the  data 
but  also  the  best  set  of  data  points. 

With  the  ability  to  evaluate  and  store  up  to  nine  sets  of 
fringe  data  on  a  cassette,  valuable  operator/computer  time  is 
not  lost  transferring  and  correcting  fori  run  inpul  data. 
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E.   CONCLUSION 

A  procedure  for  rapid  fringe  data  collection,  transfer 
and  storage  has  been  accomplished.   With  slight  modification, 
this  procedure  is  applicable  to  all  fortran  data  reduction 
techniques.   With  one  photograph  of  20  fringe  lines  and  20 
points,  the  use  of  a  digitizer  was  20  times  faster  than  that 
of  an  optical  comparator.   This  was  assuming  the  optical  com- 
parator was  used  in  conjunction  with  a  programmed  desk  calcu- 
lator to  calculate  fringe  lines  and  fringe  shifts.   The  amount 
of  time  in  the  transfer  of  data  to  the  fortran  program  was 
15  times  less  by  using  this  technique.   When  many  pictures 
are  taken  to  completely  determine  the  density  field  around 
an  object,  the  magnitude  of  time  saved  is  very  large. 

The  ability  to  rapidly  store  and  transfer  the  large 
amounts  of  data  obtained  from  a  finite  frozen  fringe  holo- 
gram results  in  a  considerable  time  savings. 
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Figure  1(c) 


Cylindrical  profile  mounted  between 
plexiglass  supports 


Figure  1(d).   Profile  and  supports  inserted  in  tunnel 
Notice  honeycomb  and  injection  tube. 
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Laser 


Lensless   Camera      /       y 


Shutter 


Hologram 
Plate  Holder 


Figure    1(e).       Schematic   of    reconstruction 
of    a   hologram. 
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Figure  2(a) 


Photograph  of  a  reconstructed  hologram 
of  a  cylinder  in  a  uniform  flow  of 
22  fps.   Notice  the  line  of  fringe 
shifts  over  the  profile  caused  by  the 
injection  of  a  thin  stream  of  carbon 
dioxide. 
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Figure  2(b) 


Photograph  of  a  reconstructed  hologram 
of  a  two-dimensional  wing  in  a  uniform 
flow  of  6  f ps .   Notice  the  Kutta  Condi- 
tion at  the  trailing  edge  of  the  profile 
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Figure  2(c). 


Graphical  output  of  potential  flow  program 
given  in  Appendix  A.   The  stream  functions 
are  for  a  uniform  velocity  of  10  ft/sec 
around  a  cylinder  of  1.18  in.  with  a  clock- 
wise circulation  of  75  ft2/sec. 
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Figure  2(d) 


Fringe  lines  of  Fig.  2(a)  superimposed  on 
graphical  output  of  stream  functions  around 
a  cylinder. 
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Figure  3.   A  view  of  the  Hewlett  Packard  digitizer, 

printer  and  calculator.   On  the  digitizing 
table  is  a  finite  frozen  fringe  hologram 
with  a  transparent  grid  and  the  electronic 
cursor.   To  the  rear  of  the  table  is  the 
acoustic  coupler. 
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data  to  he  exnanded  and  reduced,  see  Fig.  4(b) 


Figure  4(a) 


Flow  diagram  of  data  transfer  program.   Dashed 
lines  are  manual  operations.   Program  file  num- 
ber and  statement  number  is  in  the  right  upper 
corner  of  each  operation. 
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Figure  4(b). 


Flow  diagram  of  the  data  transfer  program. 
During  the  file  3  operation  the  calculator 
is  acting  as  a  terminal. 
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APPENDIX  A 
Listing  of  Potential  Flow  Program 
This  program  calculates  stream  functions  around  a  cylinder 
with  circulation. 


8U  1 
4  y  '~ 

5@  >■ 
60  't 
78  F 
80  F 
90  F 

100 
110 

120 
130 
140 
150 
160 
170 

1  P5 

180 

190 

280 

218 

220 

238 

248 

258 
p  f,  9 

278 

•-  ~7  LT 
il  f  J 

276 

2  8  8 
298 
380 
308 
309 
3 1 8 
328 
338 
348 
358 


EM****2--T.i  FLOW**** 

IM  RE  250] 

ISP  "ENTER  XMIN?  XiW,-;.  YM IN*  YMR7' 

NPU1"  71  ?  72?  VI  ?  Y2 

U  fl  L  L  7  1  :< 7  2  ?  Y 1  ?  Y  2 

R  7  I S  0  ?  +  1 

R  7 1 S  0  ?  + 1 

EM  FOR  TITLE  TYPE "LETTER" 

EM  POSITION  WITH  DISPLAY  KEYS 

EM  "EXECUTE" IS  RETURN/LINE 

REM  "STOP" TERMINATES  LETTE 


r  L"  i::.  .u 
R  MODE 


DISP 


TYPE  TITL 


LETTER 

D I S  P  ' '  I N  0  R  Y  ( .  3 )  »  2  <  <  2  5  8  > ,  P  H I  C  0  N  '•■ 

INPUT  Y6 ?  Z? 21 

D I S  P  ' '  E  N  T  E  R  P I  i;  T  U  R  E  D I  fl  &  fl  C  "i"  V  E  L  0  ( 

INPUT  R2»U1 

U=R2*U1 

R=R2/2 

T=-4*PI*R*U 

D  ISP  "  C  R I "!'  I C  f ::i !....  C  I R  C  U L  !:::l  T  1 0 N  I S : ;  T 
WRIT  4O00 

D I S P  "EN T E R  C I R C U  L fl  T 1 0  N ( C W  I S  N E 
INPUT  T 


DISP 


F  U  R  C I R  U  i....  E  G  0 


' ;:  i  i  •  I T  79k 


9  FOP 

PLOT  71 ?05+l 

DISP  "POSITIVE  PHI  BELOW  2  AXIS" 

WRIT  4008 

DISP  "ENTER  PHI > PHI  END? PHI  I  NCR 

INPUT  PjP2jP3 

PLOT  21.  ?  0  j  +1 

I  =  1 

II  =  (72  -XI).  --2 

FOR  2=71  TO  22  STEP  I ] 

Y5=Y6 

IF  7=8  THEN  358 

GOSUB  448 

RE  I  ]  =  Y 

DISP  I 


an 


i  i 
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355  1=1+1 
366  NEXT  X 

361  1=1 

362  DISP  "PHI="P 

363  HP  IT  500 

365  FOR  X=X1  TO  X2  STEP  II 

366  PLOT  X)  fl[  I  ]s  +2 
36?  1=1+1 

363  NEXT  X 

3?0  PEN 

380  IF  P=P2  THEN  420 

390  P=P+P3 

4 1 0  G  0  T  0  3  0  0 

420  STOP 

440  N=0 

4  5  0  B :::::  S  Q  R  ( X  *  X  •!■  V  *  V  > 

4  6  O  R  =  U  *  V ( U  *  R  *  R  *  V )  ■■-'  ( X  *  X  +  Y  *  V )  -  T  *  l...  0  G  ( 

465  fl= P 

470  Q=fl-P 

480  IF  RBS(fl P)  <=  21  THEN  580 

430  IE  Q*H  >=  O  THEN  51 O 

500  Y5  =  'i'5/2 

510  IF  fl<P  THEN  550 

520  Y=Y+Y5 

530  N=Q 

540  GOTO  450 

■_<  -_i  t.i   i  I  I  -.J 

568  H~Q 

570  GOTO  450 

588  RETURN 

590  P L 0  T  - R f 0  ?  + 1 

680  FOR  X  =  -R  TO  R  STEP  8 ,,81 

610  Y  =  SQR<R*R-'X*'X) 

628  PLOT  X?YH-2 

638  NEXT  X 

648  PLOT  R?0?+2 

650  FOR  X=R  TO  -R  STEP  -0.01 

660  Y=--SQR(R*R X*X> 

670  PLOT  XjYs+2 

688  NEXT'  X 

698  PLOT  -R j 8?  +2 

7O0  PEN 

710  GOTO  270 

720  END 


:pn 
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APPENDIX  B 
Listing  of  Calculator  Data  Transfer  Programs 

B.l  Listing  of  HP  File  0.   This  program  prints  the  organiza- 
tion of  the  files  in  the  tape  cassette.   By  executing  "Run 
200"  the  stored  data  of  any  picture  may  be  printed. 


18  REM  THIS  PROGRAM  BY  DAVID  POPE  USES  R  CFlLCULATERs  PRINTEF 

I  i  R  E  M  A  N  D  D I G  I  T I Z  E  R  ( H  P  9  8  3  8  i:::i  >  H  P  9  8  6  6  A  >  H  P  9  8  6  4  !:::l )  T  0  R  E  C  0  R  D  T  H  E 

12  REM  PR  INGE  SHIFT  FROM  R  PHOTO  OF  A  RE  CONS  TRUSTED  HOLOGRRI' 

13  REM  THE  INFORMATION  IS  FED  INTO  THE  OP /CMS  SYSTEM  OF  AN 

14  REM  IBM  3S0  COMPUTER  WHICH  DOES  THE  ACTUAL  DATA  REDUCTIOF 
i  5  R  E  M  T  H  E  N  0  R  M  A !....  F  L  0  W  0  !::r  D  ft  T  i:::i  I S :  F 1 1 ...  E  0 1  F I L  E  1  ?  F I L  E  2 » 

16  REM  FILE  4j  AND  THEN  FILE  3  FOR  THE  CP/CMS  HOOKUP,, 

19  FINED  O 

20  DISP  "ENTER  TIME  DAY  MONTH  YEAR"? 
30  I NPUT  Y 1  j  Y2  ?■  Y3  ?  Y4 

40  PR  I M T I  ME "  Y 1 "  DAY ,:  Y2 "  MONTH "  Y3 "  YEAR "  Y4 

5  0  F;:'  R I N  T  "  D I R  E  C  T 1 0  N  S  f  I N  D  0  R  G  i:::l  H I Z  i:::i  T  I  0  H " 

70  PRINT  "UP  TO  SO  PLANES  HUH  SO  POINTS  PER  PLANE" 

80  PRINT  "DATA  FOR  9  PICTURES  MAY  BE  GATHERED  AND  STORED" 

90  PRINT  "FILE  0: ORGANIZATION?  READOUT 

1 0  O  P  R I N 1"  ■ "  F I L  E  1 :  E  Q  U  H  T 1 0  N  S  F'  0  R  U  P  T  0  2  0  L I N  E  S " 

I I  O  P  R I H  T  "FILE  2 :  F R I N G  E  D I  S  P  L R  C  E  M  E  N  "I"  F  R  0  M  T  H  E  L I H  E  S " 
ISO  PRINT  "FILE  3: MATE  WITH  CP/CMS  SYSTEM" 

130  PRINT  "FILE  4 : LAGRANG I  AN  INTERPOLATION 

140  PRINT  "FILE  5: COMPARES  TWO  PLANES j  REDEFINES  DATR  FILES" 

ISO  PRINT  "FILE  6:  STREAM  FUNCTIONS  FOR  2-D  POTENTIAL.  FLOW" 

181  PRINT  "FILES  7  ON  ARE  DATA  STORAGE  FILES"  ' 

183  PRINT  "  STORED  DATA  READOUT  PRESS  LOR  DO  EN:.  RUNS  OO  EN" 

190  LOAD  1 

2O0  REM*****RE ABOUT  OF  STORED  DATA***** 

220  D I M  D S [ 4 1 3  ? S I C  5  ]  ?  T S  C 6 ] > 0  S E 1 5  8 ] 

230  DISP  "ENTER  PICTURE  NUMBER"? 

235  FINED  0 

240  INPUT  P 

241  PRINT  "PICTURE  NUMBER "P 

242  FINED  5 

260  LOAD   DATA  Z2jS 

270  LOAD   DATA  (Z2+1  )  »T 

271  LOAD   DATA  (Z2+2)sQ 

273  PRINT  "INTERPOLATED  DATA  FOR  CP/CMS" 
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£.  <  H 


_  _, ,-. 


284 

■"i  i~i  L." 

298 

301 
382 
383 
384 
305 
386 
387 
388 
3  8  9 
310 
311 
370 
380 
390 

400 

410 
420 
430 
440 
450 
468 
470 
471 
472 
473 
480 

500 


PRINT    "THIS    IS    FROM    PLANE "QE 2 3 

P  R I N  T    "  P  0 1 N  T    # "  ( Q  [  1  ] + 3 ) "  S  H  0  U  L  D    B  E    9 .  9  9  8  8  8 " 

FUR    1  =  1    TO    (QC 1  1  +  3) 

PRINT    I»QC I ] 

NEXT    I 

Z=SC  1  ] 

P  R I N  T    "  N  L.I  N  B  E  R    0  F    P  L  A  N  E  S    = "  S  [  1  3 

PRINT    "NO HE    OF    OPERATION    ="SE2] 

PRINT    "NUMBER    OF    DATA    POINTS    FOR    CP/CMS    -■" 

PRINT    "NUMBER    OF    VIEWS    ="SC4] 

PRINT    "SVMMETERY    OF    FIELD    ="SE53 

P  R  INT    "REFER  E  N  C  E    D  E  N  S 1 1"  Y    (  M  G  /  C  C >    =  "  T  C  L  3 

P R I  N T    "  W A V E L E N G T H    I N    M I C R 0 N S    = "IT  2  3 

PRINT    "RMAX    (CM)    =  "TL"33 

PRINT    "FILTER    ="TC4] 

PRINT    "TH    (VIEWING    ANGLE    IN    DFG>    ="TC5] 

PRINT    "DATA    SVMMETERY    ="TCG-3 

FOR    Zl=l    TO    Z 

Z2=Z1-KP*23>-14 

LOAD      DATA    Z2»D 

p  R  l   H  T    jj  j.  1  -J  „  j-,  R  T  p    p  0  j  H  T  ,;;;.    J  (, ,   p  j_  ,;;;,  ,  ,  £   „  y   ., 

FIXED    3 

I=DC 1 3*2 

FOR    K=2    TO    I    STEP    2 

J=K+1 

P R  INT    " E L E V A T 1 0 N » " D C  !<  3  j  " S I !  I F T '--■  "DC  J  3 

NEXT    K 

NEXT    Zl 

P  R  INT   "  £  ■#'  *  *  *  *  *  *  ■£  *  "£  *  "S  S  "X-  *  *  "£'  ■*■  S  *  ':!r  'i'r  -■' 

PRINT 
PRINT- 
STOP 
END 


*  ±  ir  '^  :':■  *  *  *  ±  * 
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B.2  Listing  of  HP  File  1.   This  program  calculates  the  equa- 
tions for  up  to  20  fringe  lines  in  each  picture.   The  mean 
of  the  slope  and  the  differences  from  the  mean  of  each  slope 
is  also  printed.   Two  data  points  are  entered  through  the 
digitizer  to  calculate  each  fringe  line. 


1  U  UUi'i  c~.  J.  s  J  1  )  Jc..  5  Hoi.  c.ki  J  ?  Boi.  c.'d  ..i 

i  9  P  R I H  T  "  F  0  R  C  0  R  R  E  C  T  1 0 :  ■•!  S  V  0  U 

c  b  1.1 .:.  ■:::•  r        i.J  V-.  1 1.  H  I   in  c  L  u  K  v  c  •:::• 

38  WRIT  2860 

40  DISP  "ENTER  THE  ORGIhhPRES: 


viU':' 


Mil  ■  t::..' 


J  Nil  POINT  OH  HORI 


60  DISP  "ENTER  S 

78  WRIT  2000 

80  WRITE  <9j"0 

90  ENTER  <9?*>N?Y 

100  Ji=SlW(RTN(Y/N)) 

1 10  J2~C0S<flTN(Y/;:O  ) 

120  DEG 

1 30  PRINT  "SLOPE  OF  THE  PHOTO 

140  DISP  "#  OF  FRINGE  LINES  CM 

150  INPUT  H 

160  DISP  "NUMBER  OF  PLRNES:.  MRH 

i  7  0  I H  P  U  T  2 

ISO  PRINT  "NUMBER  OF  LINES "N 

1 ?  <  0  P  R I N  I  ■  H  U  M  B  E  R  0  F  F  L  H  N 1 1  ':  2 

C&K)      J..I  ill   ilC'L  C.V   J 

2 1 8  E ~-  O 

220  FOR  1=1  TO  N 

221  DISP  "FIRST  POINT  F8R  LINE 

222  WRIT  180O 


1  5  '. 

H?s  U 


238  WRI  IF  >'■  9  s  *') 
240  ENTER  C9>*::0-: 
260  Xl=X*-J2+Y*Ji 


0  Yl=Y*J2-X-5 


U 


280  DISP  "SECOND  POINT  LIME" 

281  WRIT  1800 


■|  ••  ill-  I  .;■.•'  i- 
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299. WRITE 
360  ENTER 
320  X2~X* 
o  3  u  Y  z  ~  Y  * 
340  ML  I  ]= 
350  E=€+M 
360  R=SQR 
370  fl[  I  3= 
330  BE  I  ]= 
390  CE  I  ]= 
400  DISP 
410  WHIT 
420  NEXT 

430  M1=E/ 

431  FIXED 

432  PRINT 
440  S=0 
450  FOR  I: 
460  PRINT 
470  S  =  S-KI 
430  NEXT 
490  D  =  SQR(S.--'N) 
500  PR  IN 
510  DISP 
530  STOP 
540  LORD 


''.  9  j  ■& ) 

( 9  f  * )  X  *  Y 
J2+Y*J1 
...I2-X*J1 

(Y2-Yl)/(X2-Xi; 
E  I  3 

(ME  I 3t2+l> 
ME I 3/R 
C-D/R 

(Yl-NE I ]*Xi)/R 
"LINE  "I" 

1000 
I 

H 
4 
"THE    MEAN    OF    ' 


IS    ENTERED" 


a    TO    N 
"  L I  \ IE "  I J  "  D  I FFERENCE    FROM    MEAN    I S "  ( ME  I  ]-N  1 ) 

IE  I  3-Nl)t2 


I 


"  S  T  R  N  D  Fl  R  D    D  E  V I  R 1"  1 0  N    0  F    S  L  0  P  E  S "  D 
S  fl  T I S  F I E  D    N I T  H    S  L  0  i:::'  E  S  ? »  .  .  C  0  N  I" .  E  X 
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B.3  Listing  of  HP  File  2.   For  each  data  point  entered  this 
program  calculates  the  fringe  shift  from  the  appropriate 
fringe  line.   The  array  of  fringe  shifts  from  each  plane  are 
stored  on  tape.   A  reference  file  is  also  filled  and  stored 
on  tape. 


c-:  p  ft  m  "7  t  ..  11.  i  ■-.:■ ..  ft  ■:::■  r  o  ft  "i .,  o  •::■  r  ■'::■  ft  "i ..  r  ■  >:::■  r  ~: 

16  DIM  GS[ 41 3 

1 5  P  R I N  T  ' '  C  0  R  R  E  C  "I"  E  Fl  C  H  P  L  R  H  E  B  E  F  C 

2  0  P  R I N  T  "IF  P  0 1 H "!"  C  0  R  R  E  C  T 1 0  H  "i"  H  E 

2  5  I)  I S  P  ' '  C  0  R  R  E  C  T 1 0  N .  .  .  S "!"  0  P  C  0  H "!"  2 

3@  WRIT  5688 

35  DISP  "ENTER  PICTURE  NUMBER" j 

4  0  I H !::::  U  T  F' 

45  PRINT  "PICTURE  NUMBER "P 

58  Zl=l 

68  Z2-:"Z1  +  (P*23  )•-:!.  4 

65  FIND  Z2 

6  6  D I S  P  "  S  T  Fl  R  T  E  N  T  E  R  I N  G  P  0 1 H  T  S " 
67  WRIT  1506 

7  0  P  R I N  T  "  D I S  P I...  Fl  C  E  M  E  H  T  I N  F:;'  L  !:::i  N  E "  Z 
71  PRINT  "POINT  »  EUrVR'^  IG-1 
75  DISP  "NUMBER  OF  POINTS  IN  PLAN 


80 

I 

■iPUT 

N 

85 

G 

■  1   1  = 

H 

9  0 

1 

•-   1 

Ci  cr 

T 

:::  "j 

2 + Y  *  J 1 


100  NRITE 

5  EN" 

125  X1=X* 

130  Y1=Y*J2-X*J1 

135  J--J+1 

140  GCJ]=Y1 

145  J=J+1 

150  G[  J]  =  fl[  I  ]*Vi+B[  I  ]*Y1.+  CC  I 

J.  -J  -J  K  -  J  -  1 

160  FINED  4 

165  REM  G<1)  IS  THE  NUMBER  0 

170  REM  G(J)  IS  THE  SHIFT;.  G  < 

175  PRINT  I j  TRB10J "Y="GCK]pT 

188  IF  1=11  THEN  191 

185  1=1+1 

1 9  8  G 8  T 0  186 


6  j 

NTS 

[N    THE    PLFi! 

I£ 

THE 

ELEVRTIOH 

f't  -\ 

GC  J  3 

! :; ;  H 


HI  FT 


4  5 


1 9 1 
195 

208 
£65 
216 
215 
228 
225 
230 


240 
245 
250 
255 
260 
265 

•"J  O  c| 

290 
295 

300 

305 
310 
315 
320 
325 
338 
335 

400 

415 
420 
425 
430 
435 
440 
445 
450 
455 
465 
4  7  0 
505 
518 
515 
528 
525 
Fi  3fl 


540 

•J  4  -J 

;558 


568 

565 


WRIT  4880 

STORE   DATA  22? G 

IF  Z1=Z  THEN  215 

Z1=Z1+1 

GOTO  60 

D I S  P  "  S  fi  T I S  F I E  D  U I T  H  T  H  E  S  E  P  0 1 H  T  S  ? " 

WAIT  3O0O 

DISP  "  CONT. . .EX" 

STOP 

REM    EACH  FILE  CONTAINS  ONE  PLANE   ■• 

GOTO  395 

R  E 11  *  *  *  *  *  C  0  R  R  E  C  T I  0  N  S  *  *  #  *  * 

DISP  "ENTER  LINE  NUMBER"; 

INPUT  II 

WRITE  <9<*> 

ENTER  <9j*)X*Y 

Xl=X#J£+Y#Ji 

Y1-Y*J2"-X*U1" 

R=2*I1 

GCR]=Y1 

R=R+1 

GCR]=fl[  II  ]*X1+BC  II  D*Y1+CC  II  ] 

P  R I N  T    "  C  0  R  R  E  C  T 1 0  N    P  L  R  H  E "  Z 1 "  L I N  E "  1 1  j  "  Y  = 

DISP    "IF    S A T I S  F I E  D    G 0 . . . C 0  N  T    1 9 5 " 

WAIT    5O00 

GOTO    250 

R E  M    s x  * *  FILL    T H E    R E F  E R  E N C E    F I L E * * * * 

DIM   SIC5]jTSC6] 

Z2=P*23-16 

FIND    Z2 

S[  1  ]=Z 

DISP  "ENTER  MODE  OF  OPERATIONS)  "i 

INPUT  SEE] 

DISP  "ENTER  NUMBER  OF  VIEWS  a)"! 

INPUT  SC4] 

DISP  "ENTER  SYMMETRY  OF  FIELDS)"? 

INPUT  SC5] 

DISP  "ENTER  #  OF  POINTS  FOR  CP/'CMS"  5 

INPUT  SC3] 

STORE   DATA  22 ?S 

D I S  P  "ENTER  R  E  F  E  R  E N  C E  D E N S I T Y ( M G / C C ) "  3 

INPUT  TCI ] 

DISP  "ENTER  WAVELENGTH  IN  MICRONS"; 

INPUT  TC2  3 

DISP  "ENTER  RMFiNCCH)"; 

INPUT  TC3] 

DISP  "ENTER  FILTER"? 

INPUT    TLA  1 

D I S  P    "  E  N  T  E  R    1"  H  ( V I E  N I N  G    fl  N  G  L  E    I N    D  E  G ) "  5 

INPUT    TC5J 

DISP    "ENTER    DATA    SYMME"!"RY<8> "  8 

INPUT    TC6] 


VI  j  "SHIFT="GCR] 
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570    Z2=Z2 

H 

575    STORE 

DRTfl   Z2?T 

580    PRINT 

"IF    YOU    WANT    TO    REDUCE    THIS 

590    PRINT 

"IF    YOU    WANT    TO    GATHER    NO RE 

600    PRINT 

ii       n 

610    MSP 

'DATA    STORED    FOR    PICTURE "P 

620    STOP 

630    END 

:'ICTURE    THE 
DATA    FOR    Nt 


N, 


LOHD    5 
:t    PICTURE 
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B.4  Listing  of  HP  File  3.   This  program  recalls  the  array  of 
expanded  data  from  the  tape  cassette  and  then  the  calculator, 
acting  as  a  terminal,  automatically  sends  the  data  to  FILE 
FT03F001  in  the  IBM  computer. 


5  PRINT  "FOR  DIRECT  IONS 

1  b   Kc.il   DC  I   H L- U U ■-■  I  i 

28  REM  DIRE  £611  0 


It       ID    L.  I  ::;T1  y  '   I  I  8  «   ::   »  E)\" 

,-■  i~i  i  i  n  i  i""  1-4   r. :  i  i"i  I  i t  i-i   -i  .■ ."   r.  1. 1  r.   t  i  i  ft  i. 1   i "r  1. 1 


r  i...  i:::.  P. 

;.(  PRONE:  AND  PLACE:  IN  C0URI.J 

ibhK  h  hlbri  hi  I  l,  rib.  u     I  Urn:::. 

ij  THEN  "'"HE  LUURLER  IS  OH 


38  REM  RED  LIGHT  0 

3 1  R  E M  N  i-i  I  T  E  L  I  G  H  T  0  H  U  l-l  E  H  I" R  A H  S  M  I  T  T  I  H  G  /  R  E  C  1 1 

48  REM  TYPE  IN  "TERM. . .EX" 


Af 


:R.  ■■-!.:  Mo 


REI 


OCEDURE 


5G  REM 

68  REM  TYRE  LOGIN  PF 

i   d  ken  I  i  r  h    J.  r  L  lIic 

80  REM  TYPE  "TTY  FILE  FT83F80I  78" 

81  REM  LAST  #  IN  FT83F. .  SHOULD  BE  PICTURE  # 

82  REM  SHOULD  D  ISP"' READY  PAPER  TARE  IN  RE  ODE? 


REM  TYRE 


"  I"":  I"  L  I  I--  h.H  T  T 

p.-nNon  i  I 


m 


:L  CARDS  NOR' 


r  "i"  i  r 
J.  I  L...C 


98  REM  I'HIS  IS  THE  REED  IN  DATA  FILE  ROR  THE  C 

91  REM  TYRE  IN  TITLE (#1)  AND  TRANSMIT 

92  REM  TYPE  IN  TITLE  (#2)  FIND  TRANSMIT 

1 8  8  R  E  M  P  R  E  S  S  "  R  U  N .  .  .  E  X "  t  T  H I S  S  T  fl  R  T  S  T  H  E  C  A  L  i 

1 1  8  R  E  M  .  T  R  A  H  S  M I T  A  H  Y  T  H  I N  G  ' '  .0  0 1  ■!  T  R  0  L  <  F  6 )  -  S  -  T  R  Fl ! 

128    DIM    RIC5]jBSC6]jSSC  150] 

125    FINED    5 

1 3  8    D I S  P    "  E  N  T  E  R    P I C  T  U  R  E    N  U 11 B  E  R "  5 

140    INPUT    P 

150   Z2=P*23-16 

160    LOAD      DATA    22? A 

170    LOAD      DATA    (22+1 ') j B 

180    LOAD      DATFl    (Z2+2>?S 

183    R[3]=INKS[  1  ]> 

i  3  4    P  R I N  "I"    "  V  0  U    A  R  E    RED  U  C I N  G "  A  C  3  ] "  D  i::;l  T  A    P  0  I N  T  S 

135    PRINT 

198.  DISP    "ARE    Y00    READY    TO    SEND    DATA?" 

288    WAIT    3888 


:ULAT 
I  SUIT 


F9> 


'/CMb 
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285 
210 
226 
230 

•"  •"[  LT 

<z.  O  -J 

240 
250 
260 
265 
266 
270 
280 
290 
295 

300 

310 

O  J.  J 

320 
380 
348 
358 
368 
388 
390 
395 

400 

401    ; 

402 

403 

410 

420 


fl[3]= 

DISP 
STOP 
NRITE 

PRINT 

wait 

WRIT 

WRITE 

PRINT 

PRINT 

WRIT 

WAIT 

G=INT 

L='G+i 

FOR    I 

I  l  r~i  T  T  l~ 
W  K  i    I   C 

PRINT 

WRIT 

WRIT 

FORMA 

FORMA 

FOR  MR 

NEXT 

WRITE 

FORMR 

PRINT 

PRINT 

PRINT 

PRINT 

STOP 

END 


IHTC3C  1  ]) 

"IF    SO    j  PRESS    CGI 


.  b! 


( 4  j  340;.'  FIT.  1  li  ,  HE  2  1  j  AC  3  ] >  RC  4  3 ,  RC  5  ]  >  :l  9 

fl[  1  35RE235AE335AE435AE53 
3888 
3888 

( 4 »  350 ) BE  1  ]  ,  BE  2  ]  j  BE  3  3 1  BE  4  ] »  BE  5  3  j  BE  t 

BE  1  3 j  BE  2  3?  BE  3  3  5  BE  4  3s  BE  5  J?  BE  6  3 

3000 
3000 

( (AC  3  J.-'6)-!l  ) 


19 


=1  TO  (0*6)  STEP  6 
(4j368>SE  I]jSC  It! 
SE  I  3  5  SE  1  +  1  35  SE  1+2 

3000 
3000 

T  5F8.0jB 

T  6F10.5?B 

T  SFH.SjB 

I 

'■.  *r  ?  ■-'  ?  -J  .■'  J.  ~'  3   J.  _■  !'  1  :::' "  J.  -:'  !'   .i.  "■  ?   .!.  1/ 

T  B  j  B  s  B  j  B  j  B  ?  B 
"DATA  SENT  TO  FILE  IN  CP/CNS" 
"TO  REDUCE  DATA (EXECUTE  POPE! 
"OUTPUT  IN  FILE  FT88F881" 


5  SE  1+3  3?  SE  1+4  35  SE  1+5  3 


*    POPE! 
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B.5  Listing  of  HP  File  4.   This  program  expands  the  fringe 
shift  data  by  using  Lagrangian  interpolation.   The  array  of 
expanded  fringe  data  is  stored  on  the  tape  cassette. 


IS 

15 

16 
•t  -■? 
j.  i 

2  8 

30 
40 
50 
60 

70 
80 
90 

100 

i.  1 0 
120 
190 

200 

210 
220 
230 
240 
258 
288 
278 
288 
298 
388 
328 
338 
348 
358 
368 
378 
374 

•". "?  cr 
O  i  -..i 

376 
388 
398 
488 
4  1  8 
415 
428 
438 
448 


in 

m  y,i 

i[20]j 

YSC 

2@]jD: 

:;[  4 

PR 

INT 

"'10   £ 

TOP 

PR  IN" 

"IN! 

PR: 

INT 

Wfl 

I  1    I 

iyuy 

ni 

SP    ' 

ENTER 

PI 

JTURE 

!:!•:     i 

IN 

PUT 

PIjZS 

ni 

SP    ' 

ENTER 

# 

OF     PI! 

'.■  1  Ui' 

IN 

PUT- 

Q 

PR 

INT 

"#    OF 

PI 

JTURE:: 

"1  j 

H= 

U ~:  8 

p= 

R  =  UI 

i  1=1 

N= 

37=0 

r 

G 

OR   i 
[  I  ]= 

=  1    TO 
:  9 .  9  9 

1  ■..'• 

-j 

1 1 
n 

EXT 

i" 

f::' 

RIN" 

"    "PIC 

TUR 

E "  P  1 ,  : 

PIT 

7 

-™i "J  •' 

:+Pl*2 

"J  —  1 

.-j 

L 

8  fill 

jjflTfl 

c'L  o 

?  D 

1/ 

r-. 

'l 

=DC  : 

] 

1  ]»GS[ 


150 

.  DEL 


J  U  I  C 

5  2  8  j 


:8] 


J  d  ::::.  K  fc.  J..'  U  L  I::. .:..!  J.  o  U 


i  -  o 


X=DC  J  2  ] 
Y=DC J2+1  ] 


•, 


H=FNZ1 

IF  N<K  THEN  258 

IF  H#8  THEN  415 

GOSUB  6GQ 

DISP  "#  OF  POINTS 

INPUT  M 

FOR  1=1  TO  M 

DISP  "USE  POINT  #" 

INPUT  UC I ] 

FINED  8 

PRINT  "POINT" LIE  I  ]" 

FINED  4 

NEXT  I 

U=2 

DISP  "MIN  ELEV. j MR 

INPUT  fl?»B7»S8 

PRINT  "INTERPOLRTI 

PRINT  "MIN  ELEVRTI 

PRINT  "URN  ELEVRTI 

PRINT  "«  OF  POINTS 


J    USE    ( 


SED' 


rj 
0H"f 

0N"E 


3RSEI 

17 


II    Of 
ON' 


POIN 
"DAT 


rsM; 


■■I  i  •:.:< 


460  GC2  3  =  Z2 

470  11=2 

488  S7=<;B7-fl7 >.■•■':: 88-1 

.'i 

588  FOR  X=R7  TO  B7  S" 

IEP 

510  GC 1 1  +  1  ]=FNRX 

515  rusp  <I1-1) 

520  PRINT  (Il-D-^GI 

:  in 

525  GC  1  3=  I  1  —  1 

538  11=11+1 

548  NEXT  X 

598  STORE   DAT fi  (PI* 

22~] 

591  PRINT  "***##"GE 1 

J  "PC 

592  PRINT 

600  H  =  H+1 

610  IF  H=Q  THEN  640 

628  P1=P1+1 

630  GOTO  88 

£40    i.:isp  ■' END  OF  LH. 

I H 1 

641  WRIT  4088 

642  PRINT  "TO  REDUCE 

Dfil 

643  PRINT 

658  SI  UP 

,.-  .-  .-•.        r"i  HTM  "i"    ,!  D  i'"i  T  L  !  T   JJ.  " 

o  b  U  r  r-.  i  h  i    r  u  J.  h  i  it  ? 

665  FIXED  4 

670  FOR  1=1  TO  H 

680  Print  i  >  xc  i  d?  yc  i 

"i 

690  NEXT  1 

695  PR  IN! 

708  RETURN 

718  END 

?  L-:'  8  D  L 1  -  r  H  Z '-.  X  ) 

730  FUK  I=N  !U  1  y I  EI 

"'  --1 

748  IF  X<XL  1  ]  THEN  7' 

?8 

750  XL  1  +  1  ]=X 

I  O  O   1  L.  J.  +  X  -1  —  1 

778  RETURN  N+l 

780  XL"  i  +  i  :i  =  xc  I  ] 

790  YE  1  +  1  ]=YC  I  ] 

880  NEXT  I 

8 10   1=0 

820  GOTO  750 

830  END 

840  DEF  FNFKZ) 

868  S=0 

988  FOR  12=1  TO  N 

998  T=l 

1008  I=UC 12  3 

1010  FOR  J  1  =  1  TO  M 

102O  ,J  =  IJ[  ,11  3 

1  O  3  O  G 0  S  U B  1080 

1  3 


IN" 


filch    HMD 


.OHJJ 
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i  0  4  8 

NEXT  Jl 

1 6  5  0 

S=S+T*YC 

I  ] 

1060 

NEXT  12 

1070 

GOTO  112 

;0 

1030 

IF  XL  I  3= 

=XC  J  3 

THEN 

111 

3 

1090 

T=T*-<Z-> 

:c  J3>- 

'<XC  I 

]-XE 

j  3: 

1110 

RETURN 

1115 

T=B 

1120 

RETURN  £ 

1130 

END 
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B.6  Listing  of  HP  File  5.   The  beginning  of  this  program 
compares  two  pictures  at  one  plane.   The  second  section  of 
this  program  redefines  all  of  the  data  files. 


16    RE M    C 0 M P A R E    2    P L fl H E S ,  .  „  R U N ?  E N 

2  0    R  E  M    R  E  D  E  F  I  H  E    I)  fl  T  Ft    F I L  E  8 .  .  .  R  U  H  3  4  9 »  E  X 

3 0    R  E  M    P  0  T  E M  T !  fl  L    F  L.  0  N    0  V  E  R  ' C  V !....  I: N  D  E  R .  .  .  R  U  N  5  ?  i 

40    DISP    "COM PARE    THE    TWO    PLANES" 

50    DIM   ESC  41  LFSE41  ] 

bO    WAIT    300O 

70  DISP  "ENTER  THE  PLANE  NUMBER"? 

80  INPUT  Z 

90  DISP  "ENTER  PICTURE  NUMBERS"! 

100  INPUT  P1,P2 

1 1 0  P  R I N  T  "  P  L  fl  HE"  Z  "PI  C  T  U  R  E  S "  P  i »  P  2 

128  Zi=Z+<Pl*23>-l4 

130  Z2=Z+(P2*23) 14 

140  LORD   DATA  Z'inD 

150  LOAD   DATA  Z2?E 

ISO  FINED  4 

178  IF  DC  1  KEC  1  ]  THEN  230 

ISO  R=EC 1  1*2+1 

ISO  FOR  1=2  rO  R 

2O0  FL I ]=DC I  ]-EC  I  ] 

210  NENT  I 

220  GOTO  £70 

230  R=D[  1  ]*£h  1 

240  FOR  I  ==2  TO  R 

250  FE I ]=DC I 3-EC  I  1 

260  NENT  I 

270  FOR  1=2  TO  R  STEP  2 

280  -J  =1  +  1 

2  9  0  P  R I H  T  "  D I F  F "  F  C  I  3 »  T  Fl  B  2  0  j  F  E  J  ] 

30O  NENT  I 

310  STOP 

320  GOTO  40 

330  END 


!,•■; 


340 

REM  i 

J*** RE DEFINES 

350 

FIND 

I'' 

360 

FIND 

r 

370 

FOR  F 

:'=1  TO  10 

330 

Z2  =  p4 

=■23-16 

390 

FIND 

i^l  £Z. 

400 

LISP 

"PICTURE  IS  : 

410 

WHIT 

■j08 

420 

MARK 

1)5 

430 

FIND 

( Z  2  +  i  ■' 

440 

DISP 

"PICTURE    "F 

450 

WRIT 

208 

460 

m  n  n  !.-■ 

1j  12 

470 

FIND 

■-  i-'—  £-.    '    i.*.-  .-■ 

430 

DISP 

riLi  1...1  r:.  b 

430 

WAIT 

288 

500 

MARK 

1  ?  388 

518 

FIND 

(22+3) 

530 

DISP 

PICTURE'T 

538 

MARK 

20?  32 

540 

NEXT 

P 

550 

STOP 

560 

END 

DATA    FILES**** 


54 


APPENDIX  C 
Printed  Results  of  the  Calculator  Programs 

C.l  Results  from  File  0,  File  1,  and  File  2. 


TIME 
DIREC 

UP  TO 
DATA 
FILE 
FILE 
FILE 
FILE 
FILE 
FILE 
FILE 
FILES 
ST  OR 
FOR  C 
SLOPE 
HUM  BE 
HUM  BE 
THE  M 
LINE 
LINE 
LINE 
LINE 
LINE 
STAND 
CORRE 
IF  PO 
PICTU 
DISPL 
POINT 

1.  80 
2.88 

3 .  8  8 

4.  88 
•j  ■  y  ►j 

DISPL 
POINT 

1.80 

2.  00 

3 .  O  0 

4.  0O 


1110 

TIONS 
28  PI 
r  u  K  y 
8:  ORG 
1 :  EQU 
2;FRI 


4?  LHL 

■J'  Ln ij 

6"    CO 

~'         i~   L  I 

i        LM-i 

ED    DO' 
ORREC 


D 

AND 
FINE 
P 
N 
T I  0 


1 1 


1 .  088 

2  „  0O0 
8.  000 

4.  800 

5.  000 

RRD    D 

CT  EH 
INT  C 
RE  NU 
RCEME 

tt 
yO 

00 
00 
00 
00 
HOE  ME 

tt 

00 
80 
0  8 
00 


AY  14    MONTH  5     YEAR  1974 

ORGANIZATION 
3  WITH  28  POINTS  PER  PLANE 
TURES  MAY  BE  GATHERED  AND  STO 
AT  I  ON?  READOUT?  COMPARL  2  PIC 
NS  FOR  UP  TO  20  LIMES 
DISPLACEMENT  FROM  THE  LINES 
WITH  CP/'CMS  SYSTEM 
IAN  I N  T  E  R  P  0  L  A  T I  0  N  ( N  0  K  M  A  I...  L  Y  N  0 
IAN  INTERPOLATION  USED 
E  TWO  PICTURES  AT  A  PLANE 

DATA  STORAGE  FILES 
ERDOUT  P P E S - :. i  L  0  A D 8  E  ■%  ?  R U r -1 2 & 8 
S  YOU  MUST  START  OVER 


:ED 
UR 


i  J  '6  :  U  '■' 


f: 


OF   HE  P H 0 T 0  i  S  9 . 7 6 8 3 9 1 
R  0 
R  0 

LAN 


DEGREES 


E  SLOPES  IS    --TL0315 
IFFERENCE  FROM  MEAN  IS  8.0041 

I F  F  E  R  E  H  C  E  F  R  0  M  M  E  Fl  N  I S 0.00  9  5 

IFFERENCE  FROM  MEAN  IS  0.0273 
I F  F  E  R  E  N  C  E  F  R  0  M  M  E  A  N  I S  -  0 .  O  3  3  0 
IFFERENCE  FROM  MEAN  IS  0.8111 
HON  OF  SL  OPES,  9.  8283 
LANE  BEFORE  MOVING  TO  NEXT  PL. 
CTION  THEN. . .STOP  CO NT  258 


IT  I 


ELEVATION 

SH 

Y=  8.3156 

8. 0O85 

Y=  8. 3567 

0  .  0  o  i  6 

Y=  8. 3863 

8. 0378 

Y=  8.4175 

0. 0330 

Y=  8.4372 

8. 0722 

N  PLANE  2.0O0O 

ELEVATION 

SH 

Y  -  0 .  2  5  9  7 

8. 0537 

Y=  0.2893 

0. 8528 

Y=  0.2992 

8. 1233 

Y=  O.3107 

0. 1  146 

FINE 


IF 


I  FT 
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DI'S PLACEMENT  IN  PLANE  3.  6060 

POINT  tt  ELEVATION 

1.000O  Y~  0.232? 

2.00OO  Y=  6.3238 

3. 0006  Y=  O. 3650 

4.0060  Y=  O.4044 

5.80O0  Y=  O. 4258 


SHIF" 

O.0198 
6. 8022 
O. 0559 
■-  0 .  8  8  4  7 
0.0677 


I F  Y  0  U  W  fi  N  T  T  0  R  E  D  U  C  EI  T  H I S  P I C  "i"  U  R  E  "!"  H  E  N .  ,  .  L  0  R  D  5 
IF  YOU  WANT  TO  GATHER  MORE  DATA  FOR  NEXT  PICTURE, 


l.i  o  i  ur  r  r-.  i  N  i  1  Hla  .  .  .  .  jjfc.L 

:  OF  PICTURES  TO  BE  REDU 

'ICTURE  3  PLANE  2 

'01  NT  #  X 

1 .  y  y  y  y  0 .  z  5  y  / 

2. 0088  8. 2893 

3.880O  0.2992 

4. 8888    •  6. 318? 


i  Z  8 : 


i 


0,,  852o 
O.,  1233 
8. 1 146 


POINT    1 

USED 

POINT    2 

USE 

~  rt 

POINT    3 

US! 

"ID 

POINT    4 

MfiF 

rn 

IN  I'FRPOLH' 

"  1 0  N 

B 

R  S  E  D    0  N 

K-i  t  l  i      r-  1    r-  i  i  ri " 

II  J.  ll      L".  LLVH 

"ION 

8 

.  2588 

MAX     EL  EVA" 

"  1 0  N 

8 

.  3268 

it    L.i  F"     P  U  1  r-i 

O      1  • 

'.'     B 

A  i-i  tg$  ij 

i .  yyyy 

8. 2566 

2. 0888 

6.2558 

3 . 8  y 8 y 

6.2617 

4. 8868 

t',1  c  c.  t<  f  -J 

5. 8608 

B-27 3 3 

6. 8888 

I.1    B      L_     t          ,.■     [.... 

7, yyyy 

6. 2858 

i:::u  0660 

6. 2908 

9. 8868 

8.2967 

18. 8668 

8. 3825 

1  1  .  0008 

8.38 1:'  3 

12. 8888 

8.3142 

13. 6688 

6. 3268 

-;-  +•  -(■  +•  -i-       i    -|      ( 

1 8  0  6 

POINT 

i  6  0 !  7 


DATA  PU If 


0.3274 

8. 1389 

8„ 6223 

-8. 6354 

8„ 6478 

6  ,,6254 

8 ,,  6 1 6  4 
0. 3657 
0. 1894 
6. 1348 
6. 1289 
0  „  8  7'  8  9 

0. 0282 

GENERATED  AND 


•IT 


TOI 


0   REDUCE    DATA 
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C.2  Printed  output  of  readout  of  stored  data.   Listing  of 
stored  data  obtained  by  executing  "Run  200"  on  File  0. 

PICTURE  NUMBER  3 

I N  T  E  R  P  0  L  fl  T  E  D  D  fi  T  R  F  0  R  C  P  /  C  M  S 

THIS  IS  FROM  PL RUE  2.08060 

f :'  0  INT  #  16.60000    S  H  0  U  L  B  B  E  9  9  9  9  9  9 .  0  8  8  0  0  8 

1  .  00880  3.3.  00800 

2.00880  2.00000 

3.08888  8.32742 

4.  00000  0. 13893 

5.08880  8.82238 

6.  88888  -6. 03539 

7. 88888  -8.  84788 

8.80088  -8.02544 

9 . 8  0 0  8  8  8.81642 

18.00000  0.66567 

1 1 .  00886  6. 18943 

12. 88886  6. 13431 

13. 88686  6. 12893 

14.8688  6  8 . 8 


-i  I   c 


15.86886  -8.82821 

1 6 .  6  8  8  6  8  9  9  9  9  9  9 .  8  8  0  6  8 

L  t  i  I  m  ri  :-■  i",      fi  r-      n  I    i'"r  u  i  i~  r-      __      •-       i-:.  i-:t .-.  i":t  t":i 

i  -i  U  FI  b  c  i":.    U  i "     r  L  h  1 1  c  o    -    o  .  0  d  u  k>  U 

M  0 1)  E    0  F    0  P  E  R  fl  T 1 0  N    =    3.8  8  8  6  0 

N U  M B  E  R    8 F    D R  T  R    P 8 1 N  T  8    f:  0  R    C  P  ■••"  C  M  S    =    1 3 .  6  8  8  8  8 

NUMBER  OF  VI  ENS  =  1 . 86868 

S  V  M  M  E  T  E  R  V  6  F  F I E  L  D  =  8 .  6  8  8  8  6 

R  E  F  E  R  E  N  C  E  D  E  H  S I T  V  ( M  G  /  C  C )  =  1 .  7  4  5  8  8 

N  fl  V  E  L  E  N  G  T  H  I N  M I C  R  0  H  S  =  8  „  6  9  4  4  8 

RMfiX  (8M)  =  2.78O00 

FILTER  *  0.00000 

I  H  ( V I E  N I H  G  R  N  G  L  E  I H  D  E  G )  ::::  8 .  6  8  8  8  6 

DflT.fi  SVMMETERY  *  8.88888 

5.88686  D  fl  T  fl    P  6 1  H  "I"  S    I N    P  L  fl  N  E    1 .  8  6  8  8  6 

E  L  E  V  R  T 1 6  N  *    8 .  3 1 6                                   S  H I F  T  =  8 .  6  6  3 

E  L  E  V  fl  T 1 0  N  *    0 .  3  5  7                                   S  H  5:  F  T  =  8 .  8  8  2 

E  L  E  V  fl  T 1 0  N  =    0 .  3  8  6                                   S  H I F  T  *  8 .  8  3  8 

E L  E V fl  T 1 0  H =    0.41 3                                  8 H I F T *  6.8  3  3 

E  L  E  V  fl  T 1 0  H  =    0 .  4  3  7                                   S  H I F  T  *  0.0  7  2 

4.O0  0    D  fl  T  Fl  P  0 1 N  I  S  I H  P  L  fl  N  E  2 .  8  8  8 

E LEV AT  I  OH =  8.268              SHIFT*  6.654 

ELEVATION-  8.239              SHIFT*  6.653 

ELEVATION*  8.299              SHIFT*  8. 123 

ELEVATION-  6.311              SHIFT-  6. 115 

5.86  6    D  fl  T  fl  P  0 1 N  I  S  I N  F'  L  fl  N  E  3 .  0  8  8 

ELEVATION*  8.283             SHIFT*  8.820 

ELEVATION*  8.324              SHIFT*  8.862 

ELEVATION*  8.365              SHIFT-  6.656 
ELEVA1 ION*  6.484             SHIFT*-0. 085 

I  I  I  VA1 ION*  6.426             SHIFT*  0.068 
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APPENDIX  D 
Fortran  Subroutine  Which  Inputs  Data 

This  subroutine  is  in  the  Fortran  data  reduction  program, 
The  title,  reference  numbers,  and  fringe  shift  data  are  read 
from  FILE  FT03F001. 


SUBROUTINE  DATAC I T, QA, QB, QS) 

******  THIS  INPUTS  DATA  FROM  FILE  FT03F???  ****** 


REAL  IT(12),QB(6),QS(150) 
I NTEGER  QA(5) 
READ(3,981)( IT( I ), 1=1,12) 
WRITE(6,10G0)( IT( I ), 1=1,12) 
READ(3,982)(OA(l ), 1=1,5) 
WRITE(6,1001)(QA( I ), 1=1,5) 
READ(3,983)(QB( I ), 1=1,6) 
WRITE(6,1002)(QB( I ), 1=1,6) 
ID=6*(QA(3)/6+l) 
DO  980  1=1. ID. 6 
READ(3,98U)(QS(  l-f-J-1),  J=l,6) 
WRITE(6,1003)(OS(I+J-1),J=1,6) 

100  0  F0RMAT(1X,6A8) 

1001  F0RMAT(1X,5I8) 

10D2  FORMAT(1X,6F10.5) 

1003  F0RMAT(1X,6F11.5) 

980  CONTINUE 

981  FORMAT(6A8) 

982  FORMATC5I8) 

983  FORMATC6F10.5) 
9  84  FORMAT(6F11.5) 

IX=QA(3) 
DO  985  I =1, IX 
QS( I )=QS(l+2)  . 
98  5  CONTINUE 
RETURN 
END 
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